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Abstract 

We report the first attempt to mode decomposition of the velocity field 
using wavelet transforms in the application to the magnetized turbulence 
analysis. We compare results of the separation with respect to the global 
and local mean magnetic field for models with a strong and weak external 
magnetic fields. We show, that for models with a strong external magnetic 
field our results are relatively consistent for both methods. However, 
for superAlfvenic turbulence with spatially variable direction of the local 
mean magnetic field the differences are substantial. 



Wavelet transform 1 of a one dimensional (ID) signal s(x) produces a 
function S(a, b) of a dilatation or contraction parameter a and a translation b 
calculated from the equation: 



s{x)^S(a,b) = a- 1/2 J i/j f^-H s{x)dx, (1) 

where ip is an analyzing function called a wavelet. The wavelet transform can 
be also performed in the Fourier domain 



S(a,b)=a 1/2 ^{ak)s{k)e ibk dk. (2) 

J — oo 

In general, the function if) can be arbitrary. However, for usefulness of the 
transform, we use functions of particular properties. The function ip is assumed 
to be well localized both in space and frequency domains. In addition, tp must 
satisfy the admissibility condition: 

|V3(fc)| 2 §<co. (3) 
> \ K \ 
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This condition guarantees the reversibility of the wavelet transform. In most 
cases this can be reduced to the requirement that if) has zero mean: 



c,/, = / %jj(x)dx = 0. (4) 

J — oo 

In addition, we can require for tp to have a certain number of vanishing moments: 

x n ij){x)dx = 0, n = 0, 1, N, (5) 



which improves the efficiency of ip at detecting singularities in the signal. 

The transformation s(x) \— * S(a,b) may be inverted exactly, which yields to 
a reconstruction formula: 

•(*) = jT db£ ^ S(a, b). (6) 

The mode decomposition is usually considered with relation to simulation 
data, which in the case of turbulence are more than one dimensional. To apply 
the wavelet transform, we have to define the transform in the multidimensional 
form. Following [3] we can define multidimensional wavelet transform as 

S(a, b) = a- n ' 2 J V {^-^) s ( x ) d " x > (7) 

where x and b are n-dimcnsional vectors of position and translation, respec- 
tively. We assume here an isotropic wavelet ip(x) — ^'(|x|). The transform in 
Fourier domain and inverse transform can be expressed in analogous way from 
Eqs. ([2D and ©. 

The function S(a,h) can be interpreted here as an amplitude or energy 
contained by a wave packet of dilated or contracted by parameter a at position 
b. This interpretation allows to treat a vector field like a set of three scalar 
fields, so its wavelet transform is simplified to three transforms of each field 
component. 

We present results obtained for two 3D numerical experiments of compress- 
ible magnetohydrodynamic (MHD) turbulence for the Alfvenic Mach numbers 
Ma about 0.7 (Bo — 1.0) and ~ 7 (Bo = 0.1). Bo is the initial external mag- 
netic field. Both models have the same initial pressure which corresponds to 
the sonic Mach number equal to about 0.7. We understand the Mach number 
defined as a mean value of the ratio of the absolute value of the local velocity V 
to the local value of the characteristic speed a or Va (for the sonic and Alfvenic 
Mach numbers, respectively). Both models are calculated with resolution 256 3 . 

We start from the model with a strong magnetic field B = 1. We expect 
that both decompositions, with respect to the global and local mean magnetic 
field, should give approximately the same results. The Alfven mode is perpen- 
dicular to the mean magnetic field, so its X-component should be significantly 
smaller, because Bo has only X-component different from zero. We present cor- 
responding vectors of the field at three intersections (XY, XZ and YZ) in Figure 
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[TJ We notice, that our expectations were correct. In the top and middle plots 
on the left column we see, that X-component of the field is zero everywhere. 
These plots are for the decomposition with respect to the direction of the mean 
magnetic field taken globally. What should we expect if we take into account 
the direction of the local mean magnetic field? In this case, we expect results 
very similar. Indeed, looking in the top and middle plots on the middle column 
of Figure Q] we see that the X-component of the Alfven mode is marginal. To- 
tal energies of the particular components for decomposition with the respect to 
the global mean magnetic field are: 0.0, 0.06, 0.07 for X, Y and Z-component, 
respectively. For decomposition with respect to the local mean magnetic field 
these energies are: 0.002, 0.1, 0.09. We see that in the latter case, the energy of 
X-component is about two orders of magnitude smaller then for the other two 
components. 

The agreement in the strength of individual components, however, is not 
enough to prove, that both methods are consistent. We should expect some 
similarities in the pattern of the fields. Again, we compare plots for correspond- 
ing intersections for both types of decomposition. We see, that the patterns are 
not reconstructed perfectly. However, the most of the structure is similar. If we 
include in comparison the plots from the right column of Figure [T] presenting 
the magnetic field for the same intersections, we see, that those parts of the 
fields which are similar in plots of the left and middle column correspond to the 
places where the magnetic field is strong and approximately directed along X- 
direction. For YZ-intersection, which shows a plane perpendicular to the global 
mean magnetic field, the correspondence of fields obtained by two methods is 
in very good agreement. 

We showed that in the case of turbulence simulations with a strong magnetic 
field both methods give quite similar separation of the Alfven mode. This 
signifies, that the assumption of the direction of the global mean magnetic field 
in decomposition can be applied for velocity field separation. However, could 
it be still justified in the case of turbulence with a weak mean magnetic field? 
The answer seems to be 'no'. Figure [5] shows the comparison of Alfven modes 
obtained with the help of both methods from the superAlfvenic turbulence. 
The comparison shows clearly large differences in the field patterns between the 
plots on the left and middle columns. If we include the local mean magnetic 
field in the separation procedure, we will obtain the X-component of strength 
comparable to the other two components: Y and Z. For separation with respect 
to the global mean magnetic field we get energies of X, Y and Z-components 
equal to about 0.0, 0.1, 0.1. After taking into account the local mean magnetic 
field, these energies are 0.04, 0.1, 0.08, respectively. 

Comparison of plots of the Alfven mode separated with respect to the lo- 
cal mean magnetic field to plots of magnetic field reveals a tendency, which 
should be expected. Precisely, the vectors of the Alfven part of the velocity 
field are perpendicular to the vectors of magnetic field at the same places. This 
is additional confirmation for the validity of the method based on the wavelet 
transform. 

We have presented the separation of the velocity field into the Alfven mode 
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Figure 1: Alfven mode for the subAlfvenic turbulence. We show vectors along 
three intersections of the model domain: XY, XZ and YZ (top, middle and 
bottom rows, respectively) . Left column shows the Alfven mode separated with 
the respect to the global mean magnetic field. Middle column shows exactly 
the same intersections, but for the Alfven mode separated with respect to the 
local mean magnetic field. For comparison we also show the magnetic field on 
the right column. 
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Figure 2: Alfvcn mode for the super Alfvcnic turbulence. We show vectors 
along three intersections of the model domain: XY, XZ and YZ (top, middle 
and bottom, respectively). Left column shows the Alfvcn mode separated with 
the respect to the global mean magnetic field. Middle column shows exactly 
the same intersections, but for the Alfven mode separated with respect to the 
local mean magnetic field. For comparison we also show the magnetic field on 
the right column. 



based on the wavelet transform in the application to the magnetized turbulence. 
We compared the validity of this method for models with a strong and weak 
external magnetic field. We showed that the separation with respect to the 
global mean magnetic field can be still valid when the external magnetic field is 
significantly larger the the strength of its fluctuations. However, in the case of 
violating this condition, this assumption cannot be valid anymore. In this case 
we are obligate to use methods based on wavelet transforms. 
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